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ABSTRACT 

We explore a possibility of measuring the absolute magnitude and the 
nature (Majorana vs Dirac) of neutrino masses, by using a novel process of 
neutrino pair emission from metastable excited atoms. Except lepton num- 
ber non-conserving processes, the neutrino pair {yv) emission is the unique 
process to directly distinguish the Majorana neutrino from the Dirac neu- 
trino, using the interference effect of identical fermions. The small energy 
difference between atomic levels makes it easier to measure small neutrino 
masses as indicated by neutrino oscillation experiments. The crucial point 
is how to enhance the rate of pair emission without enhancing the radiative 
decay. We discuss two particular cases; (1) laser irradiated pair emission 
from metastable atoms, and (2) microwave irradiated emission from circu- 
lar Rydberg states. A new mechanism of the parametric amplification to 
enhance the neutrino pair emission is pointed out when Rydberg atoms are 
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irradiated by microwave, while the radiative process may be inhibited by the 
cavity QED effect. A great variety of measurable neutrino parameters and a 
variety of experimental methods make this investigation attractive. 



1 Introduction 



The nature of neutrino masses, along with its precise values and their 
mixing parameters which appear in the weak interaction, is of fundamental 
importance to explore physics far beyond the standard model. In particular, 
whether the neutrino belongs to the special class of neutral particles described 
by the Majorana equation, or to the usual Dirac particle we are so familiar 
with, is a central issue of great interest. 

In the present work we propose a novel approach to answer this impor- 
tant issue, the Majorana vs the Dirac particle, and suggest new experimental 
methods to do this. Moreover, we would like to suggest a method to simulta- 
neously determine absolute values of neutrino masses in the same experiment. 

The neutrino masses indicated by recent oscillation experiments suggest, 
but do not determine the hierarchical mass pattern. One tends to take a 
view that two mass scales suggested by the atmospheric neutrino and the 
solar neutrino oscillation is close to two heavier neutrino masses, with a 
small correction from the lightest neutrino; 

m 3 ~ 50meV , m 2 ~ lOmeV , mi <C m 2 . (1) 

This is the case of the normal hierarchy. On the other hand, in the case of the 
inverted hierarchy one has the mass relation; m 3 m m 2 ~ 50meV , m\—m\ ~ 
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(lOmeV) 2 , mi <C m 2 . How small the lightest mass m\ is and how much 
the heaviest neutrino of mass m 3 is mixed in the flavor state v e are both 
important questions unanswered by neutrino oscillation experiments so far. 
It is desirable for a single, well-organized experiment to be able to address all 
these questions. Indeed, our proposal directly attempts to answer all these 
problems with an extra bonus; if the method works ideally, one may hope to 
embark on the neutrino mass spectroscopy, along with determination of the 
Majorana vs Dirac particle. 

Available energy difference between atomic levels is closest to small neu- 
trino masses indicated by neutrino oscillation. Other energy scales are much 
larger; for instance the tritium beta endpoint ~ 18.6/ceV. Among others, Ry- 
dberg states [1] of a large principal quantum number n have energy difference 
to nearest levels of order, 



This makes it urgent to seriously consider atomic experiments for the neu- 
trino mass measurement, if the rate lies within the experimental reach. The 
pair emission rate scales with G 2 F E 5 with the energy E, the constant being 
the Fermi coupling Gp ~ 10~ 5 GeV~ 2 , hence it is usually impossible to have 
a reasonable rate, unless some novel mechanism of enhancement is proposed. 

In the present work we discuss two different types of enhanced atomic 
transitions; 7 + I* —>■ I** + ViVj , where 7 is either a laser photon or a mi- 
crowave photon. The inital atomic state I* is a metastable excited state of 
a long lifetime, for instance > lsec, while the final state /** has a strong El 
rate to a lower level such as the ground state. The experimental signal would 
be a detection of transition to the level I** experimentally designed vacant 
initially. For unambiguous identification of a weak interaction process such as 
this neutrino pair emission it is desirable to measure a parity violating quan- 
tity such as the rate difference between initial different circular polarizations. 
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(2) 
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In the first case of laser irradiated pair emission, one uses a resonance effect 
for enhancement. Our second, and a more ambitious proposal is to utilize 
an inherent instability and its associated enhanced decay of neutrino pair 
emission when Rydberg states are irradiated by a strong microwave field, at 
the same time using the principle of inhibition of ordinary radiative decay in 
a microwave cavity, a cavity QED effect [2] . 

The laser irradiated pair emission, perhaps using a more conventional 
experimental technique, might be a shortest route towards establishing the 
largest mass m 3 and the distinction of Majorana and Dirac neutrino. The 
Rydberg atom may lead to a more complete neutrino spectroscopy, including 
a precision determination of smaller masses and mixing angles. 

To the best of our knowledge the neutrino pair emission from atomic 
excited states, either spontaneous or photon initiated, has not been observed 
so far, or even not discussed extensively in the literature, presumably due 
to a clear lack of interest. We wish to point out here that the atomic pair 
emission is ideal for a precision neutrino spectroscopy. 

The rest of the paper is organized as follows. 

In section 1 we describe in detail how to distinguish the Majorana and 
the Dirac neutrinos. Our approach uses the 2-component formalism for both 
cases of the Majorana and the Dirac fields. Neutrinos that participate in 
the standard weak interaction of the SU(3) x SU(2) x U(l) gauge theory 
are described by chirally projected 2-component spinors. We find it most 
unambiguous and straightforward to use the 2-component spinor both for 
neutrinos and electrons under the nuclear Coulomb field, in order to clarify 
and unambiguously identify the true nature of massive Majorana neutrino. 
This is done using a representation of 4 x 4 gamma matrices as given in 
Appendix A, which is not the most popular one in the literature, but is 
explained in many textbooks such as |3J. A comparison of the 2- and the 
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4-component approaches is also given in ref. [3]. 

The most popular approach uses the 4-component spinor ip with the Ma- 
jorana condition i\P = ip, which essentially reduces the 4-component ip to 
the 2-component spinor, if ~ ip c + ip. A great merit of the 2-component for- 
malism is that it uses independent variables alone. On the other hand, the 
4-component formalism uses redundant fields constrained by the Majorana 
condtion. We prefer to use independent components alone since the neutrino 
that appears in the usual weak interaction needs two components ip alone. 
In our 2-component approach it is made evident below that the distinction 
of the massive Majorana and the Dirac cases occurs only via the interference 
term of two identical particles present in the Majorana case. 

We then present the pair vv emission amplitude and demonstrate how 
the distinction arises in the two cases. The other place where the Majo- 
rana nature arises might be in lepton number non- conserving processes such 
as the neutrinoless bouble beta decay. But in the case of lepton number 
non-conservation there can be no proof that the Majorana neutrino is di- 
rectly involved, since there might be another source of lepton number non- 
conservation.. 

In section 2 we work out kinematical factors of the pair emission from 
excited atoms. In section 3 we discuss one of the enhanced process; laser 
irradiated pair emission. How the threshold behavior of photon irradiated 
pair emission differs in the Majorana and the Dirac neutrinos is described 
in detail. We then numerically estimate the rate for this process assuming a 
standard laser flux of commercially available frequency resolution. We shall 
further illustrate how to determine neutrino parameters including 3 mass 
values and the angle 13 [4J if the proposed experiment becomes possible. 

In section 4 an entirely new process of pair emission from circular Ryd- 
berg atoms is discussed. Both the standard multiphoton picture and more 
general effect of the parametric amplification is described. The enhancement 
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factor of the microwave irradiation is interpreted in the multiphoton picure 
as the existence of a great many varieties of paths of stimulated photon emis- 
sion that bridge between the initial and the final states of energy difference 
of order of the sum of the mass of two emitted neutrinos. The multiphoton 
process corresponds to the narrow band region of the parametric resonance. 
What is more interesting is the wide band region of the parametric amplifi- 
cation which is missing in the multiphoton picture, and one may expect an 
even larger, exponential growth of the rate. We discuss the unitarity bound 
on the pair emission rate in the wide band region. 

Our basic assumption throughout this paper is that the standard elec- 
troweak theory correctly describes the neutrino interaction, while their small 
masses, either of the Majorana or of the Dirac type, are generated at a 
much higher energy scale than the electroweak symmetry breaking. Thus, 
the weak current associated with neutrinos is always taken of the V — A 
form. We warn that introduction of V + A current may drastically change 
the result presented here. 

Throughout this paper we use the natural unit of % = 1 and c = 1, and 
a ~ 1/137 and a 2 m e /2 ~ 13.6eV. 

2 How to distinguish the Majorana neutrino 
from the Dirac neutrino 

The most commonly assumed method of observing the Majorana nature 
of neutrino masses is to discover the lepton number non-conservation, as 
typically exemplified by the neutrinoless double beta decay. But, this is 
not the unique way of detection. Another, and a more direct method is to 
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exploit the identical particle effect of Majorana particles [5], since a Majorana 
particle is identical to its anti-particle. 

In the non-relativistic regime where the distinction of Majorana vs Dirac 
particles is expected to appear, plane-wave solutions for the Majorana and 
the Dirac particles might appear different, and this difference might show 
up in many places if momentum of neutrinos is smaller than their masses. 
We need a systematic way to handle the most general cases. This is the 
purpose of the present section. In the end we shall show that there exists a 
representation demonstrating wave functions common to both the Majorana 
and the Dirac cases. All other representations should give equivalent results 
to this one. 

Since a convenient account of the Majorana field, in particular of the 
2-component formalism, is missing [3J, we explain in detail and give fun- 
damental formulas related to Majorana neutrinos, in particular an explicit 
plane-wave solution for the massive Majorana particle. Our approach is 
based on the 2-component spinor, and unlike many other works, in no place 
we adopt the 4-component description. This way we believe that a possible 
complication due to a constrained fermion field of 4-component description is 
avoided. We shall prove that the interference term of the anti-symmetrized 
wave function of identical fermions is the only source of distinction of the Ma- 
jorana and the Dirac neutrinos. A great merit of our 2-component approach 
is that this simple result is an automatic consequence of our formalism. 

Fine details are relegated to Appendix A, some of which should be useful 
in different contexts. 

2.1 Majorana equation 

Lorentz invariance allows electrically neutral particles to be described by 
two component spinor equation, as pointed out a long time ago by Majorana. 
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The Majorana equation for free neutrinos is 

(id t — iv ■ V)v? = ima2<p* , (3) 

with m the neutrino mass. The 2-component spinor (p belongs to an irre- 
ducible representation of the Lorentz group, unlike reducible representation 
of the 4-component Dirac spinor. The most salient feature of this equation 
is that it contains <p as well as its conjugate (p*, thus the lepton number is 
violated, or more properly, one cannot define the lepton number. Unless the 
lepton number is violated in other places of interaction, the rate of lepton 
number violating processes is proportional to the square of the neutrino mass, 
actually some weighted average of neutrino masses squared. 
The plane- wave solution to eq.Q is given by 

(4) 
(5) 

Consistent quantization as discussed in Appendix A leads to the normal- 
ized operator form of plane-wave solution written in terms of the helicity 
eigenstate of eigenvalue h, 




<p?,h{x) = c(p , h)e ip ' x u(p, h) + c\p , -h)e ip " 



u (p, h) = - 



E p - hp I P + hp 3 
2 \\ pE p (p + hp 3 ) \ h{p x + ip 2 



(7) 



The helicity eigenstate wave fucntion satisfys (<r ■ p/p)u(p,h) = hu(p,h). 
Quantization of the Majorana field as explained in Appendix A gives the 
interpretation of c(p, h) and c'(p, —h) = (c(p, h) y as annihilation and cre- 
ation operators of Majorana particles of momentum p and helicity h. 
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2.2 Weak interaction of neutrino 



We only consider weak interaction of neutrinos with electron, since our 
subject is the atomic weak process, hence we ignore heavier charged leptons 
and quarks. 

The Majorana neutrino field appears only in the form of the projected 2- 
component spinor, tp = (1 — 75)^/2 in all weak processes. We shall also write 
down the electron field operator decomposed into the 2-component form, 
which must be done using the same representation of 7 matrices as done for 
the neutrino. 

It is convenient to use the Fierz transformed 4- Fermi form including both 
charged current (CC) and neutral current (NC) interactions; 



-j= z/ e 7 Q (l - 7 5 )z/ e e7 Q (l - 7 5 )e 
G 



^= ]T P* 7a (l - 75K e ( 7 a (l - 4 sin 2 9 W - 75)) e , 



(8) 



where sin 2 9w ~ 0.231 experimentally. The relative sign of CC and NC terms 
becomes important later. 

A care must be taken of the effect of the nuclear Coulomb field on elec- 
trons. Fortunately, to orders of a and l/m e , the result is simple since 
(E — m e — V)/m e = 0[a 2 ] terms can be neglected. The result of non- 
relativistic limit is summarized as 



2\[2G F x (vlv^e + v\av e ■ e^aej 



2 4" 



t + 



l-4sin 2 9w(l + — a- V) 



1 \- t 



a + 4sin 2 9 W — (-iV - ffxV) 



This is rearranged to 
G F 



J-/ - _£ \ ^ AO- A£ 



(9) 



(10) 




) 



(11) 



(12) 



with a a = (l,a). 



2.3 Dirac neutrino and comparison with Majorana neu- 



The relation to the familiar 4-component Dirac equation is explained as 
follows. Using a representation of the Clifford algebra that diagonalizes 75 
(its explicit form is given in Appendix A), the Dirac equation is decomposed 
into two equations for two independent 2-spinors, if and %; 



Thus, the identification by % = i (T 2f* i n the Dirac equation gives the Majo- 
rana equation, eq.(j3]). 

Physical content of these two equations appears different; only 2 helicity 
states exist for the Majorana field, a particle being identical to its anti- 
particle, unlike distinguishable particle and anti-particle for the Dirac case. 
Whether the neutrino as observed in the V—A weak interaction belongs to the 
Majorana case or the Dirac case is the unsettled question facing fundamental 
physics. 

It is important to theoretically compare the chirality-projected Dirac field 
ipD — (1 — 75)^/2 to the Majorana field. The relevant 2-component operator 
corresponding to the momentum eigenstate oc e~* pt+ip " x i s 



trino 



(id t — id ■ V)y> = mx , (id t + io ■ V)x = rn<f ■ 



(13) 



ipD = b(p, h)e ip ' x u(p, h) + d\p, -h)e 



E p + hp 



(-ia 2 )u*(p,h). (14) 



\ E p - hp 
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Here u(p,h) is given by eq.(J7|). Anti-particle creation operator S(p,h) ap- 
pears here, which is distinct from the particle creation b*(p, h). 

Comparison of the Majorana solution, eq.flH]), ([7j) and the projected Dirac 
solution, eq. (|14p demonstrates the equivalence of the two wave functions, the 
difference being the distinction of the Majorana particle and the Dirac anti- 
particle eft creation. Their distinction appears only via the identical particle 
effect of two Majorana fermions, an extra term containing c^cj = — cj&j,. 
The unique process to distinguish the massive Majorana from the massive 
Dirac neutrino is thus the pair emission i/u, in which the anti-symmetrized 
wave function appears only for Majorana neutrinos. We shall later show that 
this distinction too disappears in the high energy limit of E p ^> m. On the 
other hand, neither a single nor a pair v v (not v v) emission can tell their 
distinction even for non-relativistic massive neutrinos, although they may be 
able to determine the absolute mass of neutrinos. In this sense the double 
beta decay of two accompanying neutrinos v e v e is useless for distinction of 
the Majorana and the Dirac neutrino. There are thus only two experimental 
ways to verify the Majorana nature of neutrinos; the other indirect method 
is to verify the lepton number non-conservation such as in the neutrinoless 
double beta decay. 

2.4 Pair emission 

The idea of using the decay of unstable elementary particles to verify 
the Majorana nature of neutrinos via identical particle effects is not new; for 
instance, see [5] and [6]. The problem of this approach is a huge disparity 
of energy scales; in both cases of the rare K- meson decay K — > n + ViVj [5] 
and the muon decay yU + — > e + i/ e u fl [6], the mass difference is much larger 
than anticipated neutrino masses, and even if events of this process are ac- 
cumulated statistically, there is no sensible way to precisely determine the 
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neutrino mass. 

On the other hand, the atomic energy difference is closer to the neutrino 
mass scale; 

A£ ni , na ~13.6 e yx(^-^) (15) 

n 2 n-y 

which reduces to 27 meV An / (n / lO) 3 for \n\ — n 2 \ ~ n 3> 1 for Rydberg 
states. 

Let us first discuss the neutrino pair v 2 Vi emission of Dirac particles. Note 
that the anti-particle notation V\ is necessary only for the Dirac case. The 
Dirac pair emission oc tf(p 2 , h 2 )d\pi , hi) is governed by e l ^ 1+ ^' x times the 
current matrix element, 

3D(Pihi,P2h 2 ) = felp j"^^ ' ^ h2 ) craia2U *(P 1 ' ^0 -( 16 ) 

The pair emission rate contains a neutrino current product; 




with a& = (1 , -<?)■ 

The neutrino pair current given above is to be multiplied by the electron 
current product. After this multiplication, the helicity summed quantity for 
the Dirac case is 

E ia> ■ j'I" = Ud + + a - liPW)' 

AiAa 2 \ h x h 2 hih 2 

hiihi 2 hi i hi 2 

+2%^ ■ ^(?)t + 2(f - f ) ■ x ^) . (18) 
hi\hi 2 hii hi 2 j 

The last two quantities in proportion to imaginary parts of the current prod- 
uct may contain CP-odd effects. 
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Let us next discuss the Majorana pair emission. The Majorana pair 
emission operator, c\p 2 , h 2 )c j; (pi , hi), gives a matrix element of two anti- 
symmetrized wave functions due to the anti-commutation of the Majorana 
field. The neutrino current for the pair emission is thus e l ^ 1+ ^' x times 



3u{V\hi ,p 2 h 2 ) = ^ + |^V (pi , hi)a a a 2 u*(p 2 , -h 2 ) 

^ fi-hlly ^ 2 ' "^K^^Pi > • ( 19 ) 
To derive the rate, one multiplies the electron current jf and takes its 
square. After a little algebra, one finds for the relevant quantity of the 
neutrino part, 

ju&y=^+hx^)(i+h 2 ^-)x 

tr (l-/ l2 ^-^W Q (l-/i 1 ^i)a /3 + (l <->2) (20) 
P2 Pi 

+ 7-1 tr (1 - h 2 )a (1 - hi K + (1 2) . (21) 

i<ohiihi 2 p 2 pi 



The last two terms of eq. (1211) are the interference term proper to identical 
fermions. 

It is evident that without the interference terms oc mim 2 /(EiE 2 ) in (I2ip . 
the Dirac and the Majorana emission rates are identical, by considering an 
extra factor 1/2 for the Majorana case, which is necessary after the phase 
space integration because the same configuration of identical particles are 
counted twice. Hereafter we devide the Majorana contribution by 2, antici- 
pating this overcouting beforehand. 

The helicity summed interference term thus becomes 

Ifl^S)' -?•(?)')• (22) 
Hence the current product in the rate for the Majorana emission is 

£ \3M ■ f? = £ \jo ■ J'T + fg| (j'o e (Jo) t - f ■ (J e ) f ) • (23) 
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The interference term is CP-even. Hence the CP violating effect is identical 
in the Dirac and the Majorana cases, for the neutrino pair emission. 

As an illustration, let us work out the current product by neglecting l/m e 
terms and taking sin 2 9 W = 1/4. We call this the leading approximation of 
l/m e expansion. The spin and the orbital part of the electron current is 
separated as 

jo = S aa >ai f (pi + p 2 )cfj , f = (sVl s )°i/(Pi + P2)c\f , (24) 
MPi + ft) = (f\e-^ + ^\t) , (25) 
cg> = U* ei U ej , 4 S) = U* ei U ej - l -5 l3 . (26) 

Furthermore, we consier the spin averaged rate such that 

\Y,jt(jt) ] = SMh^ +p 2 )i 2 | C g ) i 2 , l^my = mpi +p 2 )i 2 ii 0) r- 

ss' ss' 

The result of the spin average for the Dirac and the Majorana cases is 
^2^2^\jD-j \ = \aif(pi+p 2 )\ x I ( 1 + "BB")l c ij I +( 3 - IFiHK I 

+2(|-|)-Rej e x^|cl;)| 2 ) , (27) 

ss' h\h,2 ss' hih'2 

I /- , -M2 m l m 2/ Q , (s),2 I (0),2\ / 00 \ 

= -\aif{pi + p 2 )\ ^-^( 3 l c i/l ~ l c VI ) • ( 28 ) 

The long wavelength approximation for the neutrino makes this formula 
much simpler, allowing the replacement, 

\a if (pi + ft)| 2 = |(/|e-^ + ^»| 2 -> 1 . (29) 

This is valid if the wavelength of the neutrino X u ^> atomic size, or p u <C 
Zam e /n 2 . Since our main interest is in the region of the neutrino mass, this 
means, with m v < p v , 

Zam e Z . . 

m^< -~3.7keV — . 30 

n n z 
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This relation holds in the following discussion for laser irradiated process, 
but may not for pair emission from Rydberg atoms. The condition for the 
long wavelength approximation is significantly modified for Rydberg atoms. 



3 Kinematics of pair decay 

3.1 What can be measured from the pair decay 

A merit of the process of neutrino pair emission from excited atoms is a 
great variety of measurable quantities related to the neutrino mass parame- 
ter. The neutrino pair v(Uj can be any combination ij of mass eigenstates, 
as is clear from coexistence of the charged and the neutral current interac- 
tion. From the energy threshold position one can determine a combination 
of the neutrino mass mi + rrij, while the strength of the rate gives the mixing 
parameter in the form, \Ul^J e j\. Since any pair ij (altogether 6 channels) is 
conceivable, there is a great many combinations. This is why we phrased our 
experimental approach as the Mneutrino spectroscopy. The situation is quite 
different from the neutrinoless double beta decay in which one concentrates 
only on a combination of parameters | J2i U e iU ei mf | , which however attempts 
to discover the important issue of lepton number violation. 

Ideally, one can determine all masses mi with a bonus of experimental 
redundancy. In particular, we would like to emaphasize that this is the first 
opportunity to probe the smallest mass m 1 . Furthermore, we may explore 
a possibly very small mixing factor \U e3 \, which indicates how much the 
heaviest neutrino is mixed in the flavor v e . Both are important since other 
experimental methods may have no good handle on these quantities. 
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3.2 Phase space factor 

When the pair emission occurs between 2 levels of energy difference A, 
the rate is given by the phase space factor of 2 massive neutrinos viUy The 
rate is 



lf^(A) = 8G| / ^ 1+ £ 2 -A)VEMm,z,-j1 2 (31) 

Z SS> hlh 2 



The constant c^-'s are different, depending on whether the electron transi- 
tion involves spin- flip (F) or no flip (NF). They are in the leading l/m e 
approximation, 

~2 



°ij = i s) = K u ^ - 4 ( p ) > % = i 0) = ( NF ) » ( 32 ) 



EM?r=!> eM?i 2 = i. (33) 

The kinematical factor is defined, when the matrix element is ignored, as 

= ^ £""" <<Bi Bi(A " E,)p? " m?) ((A - E,) 2 - m 2 )]' /2 -(34) 
Near the threshold A ps m» + 

/£>(A) » -L( mjmj ) 3/2 (A - m, - m,-) 2 . (35) 

Note that the distinction of Majorana and Dirac particles appears only when 
i = j, since with i ^ j two neutrinos have different masses and are not 
identical particle. 

For illustration we shall give the rate near the threshold in the leading 
approximation of l/m e , further suppressing the correlation of electron tran- 
sition amplitude with neutrino momenta, 

r S ?(A) * 0l4 O) | 2 K^) 3/2 (A - m, - m 3 f , (36) 
Tg(A) « ^(IcS?! 2 + 3|cS; ) | 2 )Km,) 3/2 (A - mi - m,) 2 . (37) 
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4 Laser irradiated pair decay 



The weak rate scales with the available energy as E 5 . This means that for 
a small available energy the rate is very small. For the case of the neutrino 
pair emission from excited atomas, the phase space factor is 

^.3,xl0- s -.(|)=, (38) 

which should be multiplied by the matrix element squared. One clearly needs 
some enhancement mechanism even to hope for detectablity of the neutrino 
pair emission. We shall discuss in this section laser irradiated pair emission 
using a resonance effect. 

Let us first estimate very crudely how much enhancement may be ex- 
pected for resonant processes. Suppose that the neutrino pair emission from 
a metastable atom of lifetime I/7 is triggered by laser irradiation of flux 
Fq. The rate for photon absorption is gFq with a the photo- absorption cross 
section, and this irradiation is effective for the duration of lifetime. Putting 
these factors together, one might naively expect laser irradiated rate of order, 
o FqT^/j, with T uP the rate of order fl38|) . 

Another important factor is the energy resolution Au of laser. Convolu- 
tion of the laser spectral function with the Breit-Wigner resonance function 
of the natural width 7 leads to a factor of order, 

^ -1.5x10-^^-1^, (39) 



7A0; Au; leV 7 ' 

with ujq the laser central frequency. It was assumed that Au; ^> 7, which is 
usually valid. We shall take a standard laser flux of order, 

W TTiTTt ^ 

F = — « 6.2 x 10 20 cm~ 2 secf 1 « 1.6 x Kr 4 (el0 3 , (40) 

ev 

and the photo- absorption cross section of order, nm 2 . The laser beam power 
P is related to the number flux Fq by P = ujqFq, with ujq = h x the laser 
frequency. 
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This leads to a typical laser irradiated rate of order 
12 c 15 — p- Is- 1 a , E 



F - 2 x 10~ 18 s _1 



5 (10 -9 ^0 



p 



• (41) 



157T 3 Au'y 7 nm 2 leV Alu Wmm~ 2 

In the following we shall give a more concrete estimate for a particular process 
of laser irradiated pair emission. 

We consider [7] laser irradiated neutrino pair emission from metastable 
ions or atoms |/*), 



7 + /* 



+ ViVj , 



(42) 



where the final state |/**) is a short-lived excited state, as dipicted in Figure 1 
[8]. Detection of |/**), presumably via radiative decay into the ground state, 
gives a signature of this weak process. A measurement of parity violating 
quantity is highly desirable for the background rejection. 




Figure 1: Atomic level structure and laser irradiated neutrino pair emission 



The intermediate state \I n ) is chosen to lie energetically below the ini- 
tial state | J*) by an amount of the paired neutrino energy E,i + Ej with 
Ei -- ~ " 



I Pi + mf. The laser frequency is tuned to the next step of radiative 
transition from \I n ) to the final state The Breit-Wigner resonance 
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factor 

1 _ 1 

(E* -E n -Ei- Ej) 2 + 7 2 /4 ~ (eZ -E n - cu) 2 + 7 2 /4 ' (43) 

with the energy conservation E* + u = E^ + Ei + Ej, gives a large enhance- 
ment at the threshold Ei + Ej = rrii + m,j, if 

E n - m, - rrij , a; ps - + m, + rrij . (44) 

Thus the laser initiated pair emisssion has the threshold at the laser energy 
of 

iu th = E„ - E* + rrii + rrij . (45) 

The threshold location cu t h is expected to be measured with a good precision 
of the laser frequency. 

The initial state |/*) must be a metastable excited state, for instance 
0[2m 3 ] above another state \I n ). The intermediate state \I n ) can be either a 
ground state, or better, another metastable state. The width factor given by 
7 2 = 7 2 + 7 2 is a sum of the initial and the intermediate state contributions. 
We assume that both of these widths are of order lsec -1 or smaller. There are 
many candidate atoms or ions of this kind. Another important assumption 
is that we prepare depletion of the intermediate state \I n ), since the laser 
excitation of this state to |/**) is the crucial signature of the pair emission 
process. 

To lowest order of the weak interaction, the rate for the pair emission v^Vj 
of mass eigenstates when a laser of flux F(u>) is irradiated, is given by 



s; m , . 8G 2 F e 2 u\{I**\x\I n )\ 2 t M,r>, , 

r « M = FM 3[(^-A /W ) 2 + 7 V4] ^ ( " " Afi) ' (46) 

fi)/A\ _ tWfA\nJP)\2 , Q|»|2 



2^(A) = ^(A)(|c^r + 3|^r), (47) 
2#(A) - 2/£(A) = 4 2) (A)(|4°>| 2 - 3|4 S) | 2 ) , (48) 
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where is defined by fl34l) . and 

A/„ = E^ — E n , Ajj = E^ — E* . (50) 

The Majorana and the Dirac difference oc fj? \oj — Afi) x (jc^l 2 — 3|c^| 2 ) . 

When integrated over a laser spectral function of the energy resolution 
Alo > 7, 

FH 2ttF 
rf %-A) 2 + 7 2 /4~^7' (51) 



the rate is 



M " 3 7 A^ /i ^ ( ^ " Afl) ■ (52) 



The dipole strength is related to the natural width, 



4-rfT 

7r = A/|(f>| 2 , (53) 



which gives 



v m (:: n _ Ifor^Trifr 

f 6vr 2 G|-^-^(co - A /4 ) . (5.1) 



WqAu; 7 

We assume that the laser tuning is complete. In this last formula three 
factors are separated; the laser quality factor F / (uqAu), the atomic factor 
7r/ 7 , and the neutrino kinematical factor f^f. 

Near the threshold, denoting the tuned frequency lj by lu, 

f»(u> - A„) = ^\cf\\m imj fl\u - A fi - m, - ™i? , (55) 

fS( U - A fi ) = ^(Slc^l 2 + |4 0) | 2 )(m i m J ) 3/2 („ - A fi - rm - m 3 f .(56) 

To compute a reference rate let us take 13/4 for the asymptotic value of 
the factors of \cij\ 2 , which gives a basic unit of the rate, 

" . » . P eV 1^,10^ 



QAn 2 l) 2 Alj^i Wmm~ 2 lu Alu 7 
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We have in mind radiative rates of order, 1/7,. ~ Ins and I/7 ~ Is. Taking 
the energy scale at 0.3eV, about 3 times the pair mass of the heaviest neutrino 
50meV, then gives the rate of order, 

1 x 10 - 2 yi flM p ^ iQ-^io-V (58) 

U (0.3eV)s Wmm-i [ uj ] Au 7 ' [ ' 

This corresponds to 1 event/day for 10 16 target atoms. With a more 
experimental effort of improvement such as the use of the resonator for laser 
irradiation, an enhancement of ~ 10 3 — 10 4 may be expected and would much 
help for the improved event rate. 

The threshold suppression is large due to the square factor (uo — A fi — 
rrii — rrij) 2 , but the rate rises towards 

2 ^(^f\ 2 + \4 ) \ 2 )^4^. (59) 
15tt V 1 v 1 1 v 1 J luJ 2Atu v ' 

This pattern repeats for each pair ij, and finally approaches at much larger 

uo > 2m 3 , 

^T ( 4 + 1) 7^' (60) 

where 9/4 comes from the spin flip term, while 1 from the non-flip term. To 
determine neutrino massses, it is necessary to fit experimental data of the 
threshold rise up to an intermediate energy range. It is then important to 
have a large statistics data with a reasonable precision in the vicinity of the 
threshold, typically away from the threshold a few to several times the sum 
mi + rrij. 

The neutrino mass spectroscopy may proceed step by step. First, the laser 
frequency dependence, for instance oc (uj — Afi—rrii—mj) 2 near the threshold, 
may be used to determine mass parameters mj. Simultaneous with or even 
prior to m 3 determination at the threshold u 3 is 3 , distinction of the Majorana 
and the Dirac cases is presumably possible at uj — Ajj ~ 0[6m 3 ] 0.3e^. 
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Once the mass determination is done, one proceeds to determine mixing 
angles by measurement of the absolute rate. For instance, the sensitivity to 
the smallest, unknown angle 9 13 is large at the threshold of uj = A/j+m 3 +mi, 
since the relevant factor has a large coefficient; 

3\c$\ 2 + \$\ 2 ~2.9sm 2 6 13 . (61) 

For the heierarchical mass pattern, m 3 + mi ~ 50meV. The rate at this 
threshold is amaller by a factor ~ 1/32 than at the the threshold 2m 3 ~ 
lOOmeV. Although smaller in the rate, the #13 measurement may be possible. 

We conclude that for precision determination of absolute values of ,(i — 
1,2,3) and 6*13, pair emissions of (^3^3) , (^3^2) , (^3^1) near their thresholds 
are channels we recommend. 

A quantitative Majorana-Dirac distinction may be much helped by noting 
the rate difference at each threshold ij; 

r£M - rg(c ) = ^-^(Icfl 2 - 3| c y | 2 ) , (62) 

l4 0) | 2 -3|4 s) | 2 = -2|C/ ei | 2 |C/ ej | 2 , fori^j (63) 

= -2\U ei \* + 3\U ei \ 2 - I , for % = j . (64) 

For a reference, we give a complete rate formula including all pair chan- 
nels; 

2G 2 Flr 

r / acu 

1.) 

F{uj ; u , Aa>) 



T M (u ) = 2 -^f\ fduoY, 0{u - A fi -rm- mj ) 
71 u o J a 



x- 



cu-A /n ) 2 + 7 2 /4 



x 



c 



(°)|2 , Q|>)|2V, , A ( ^ ^ 

+ 3|cy I )( W - A fi ) Yij [ — — , — — 



-H I - 6 \ c ij I jmrnjiuj - I\ fi ) — , — 



,(65) 



Yijie^ej) = I' ^ dxx(l - x)J(x 2 - ef) ((1 - z) 2 - e 2 ) , (66) 
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dxJ (x 



(67) 



We used a notation of laser spectral function F(u ; uo , Au), which has a 
central frequency uq and an energy resolution Au. In the case of the Dirac 
neutrino the last term mirrijZij in (1651) is missing. 
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Figure 2: Ratio of laser irradiated pair emission rate for Majorana (M) and 
Dirac (D) cases to the massless rate. Vertical lines indicate 6 threshold loca- 
tions. 

Under the assumption of Au ^> 7, this rate becomes simplified as 

T M (uj) = ^-±-^J2 d ^- A fi- m i- m j) dE 1 I(E 1 ), (68) 

7TU IAU 7 4 . Jm,i 

I(E X ) = (|4 0) | 2 + Slcfl 2 )^ - A fi - E^{El - m?)( (u - A fi - E x f - mf ) 
+^(li 0) | 2 - Slc^nm.m.V^ 2 - (u - A fi - Etf -mj). (69) 

We plot in Figure 2 the ratio of two rates; the Majorana rate T M (u) and the 
Dirac rate T D (u) divided by the rate of massless neutrino pair emission. The 
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paremeters used for this figure are 

sin 2 12 = 0.35 , sin 2 13 = 0.032 , 

mi = l.OmeV , m 2 = 9.0meV , = 50.8me\^, 

constrained and allowed by neutrino oscillation data. With this small mass 
mi the energy region in which the Majorana rate is larger than the Dirac 
rate is restricted to a small region below mi + 777.2. When mi is larger, say 
mi > 2meV, the Majorana dominance persists up to slightly above the 2m3 
threshold. The Majorana dominance is also sensitive to a value of sin 2 # i2 



These rates presented here may be an overestimate for a dense gaseous 
target, since the energy resolution in this case is governed by a usually much 
larger collisional width, 



5 Neutrino pair emission from circular Ryd- 
berg states 

5.1 Circular Rydberg states 

Circular Rydberg states [10] are highly excited; in addition to a large 
principal quantum number n, it has the highest angular momentum, I = 
\m\ = n — 1. These states have the least overlap with the atomic core of 



19 j. 




(70) 
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charge (Z — l)e. Its size (r) pa n 2 dB with a dispersion y((Ar) 2 ) pa yn 3 /2aB, 
and the average momentum (p) pa l/(nas) with a large dispersion. Thus, the 
circular Rydberg state is almost classical, and the system is approximately 
described by the hydrogen-like Coulomb potential of charge e. A great merit 
of circular Rydberg atoms is that lifetime of radiative decay is very long, 
scaling as n 5 with the principal quantum number n, and ~ 0[lms] for n ~ 25 

mi. 

The wave function ipnim of a circular Rydberg state of a principal quantum 
number n is given by 

VW- 1± (n-l)(r) OC e ±l >-^(r 2 - 22)(«-l)/2 e -r/(««B) ? (71) 

with the magnetic quantum number m = ±(n— 1), and 05 is the Bohr radius. 
For large coordinate arguments r» Rydberg states in general have the 
radial wave function of the form, 

Rnn-l{r) OC r n-l e -r/(na B ) _ (73) 

Momentum representation of Rydberg states is also useful, as shall be 
shown in the following. Appendix B describes momentum space representa- 
tion of the wave function. 

For subsequent discussion it will become important to compute a corre- 
lation integral of the initial and the final wave functions of atomic electron, 

«i/(A) = / ^e- iK -%{x)U$) , (73) 

where A = q\ + g*2 is the sum of emitted neutrino momenta. We assume that 
both the initial and the final electron states have definite azimuthal angular 
momentum components, rrii ,mf along the microwave propagation, taken as 
z— axis. Using the expansion formula, 

e -iA.£ = J2i\2l + l)^(Ar)P,(cos0') , (74) 
i 

cos 9' = cos #a cos 9 + sin 9& sin 9 cos(<^a — V 9 ) > (75) 
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along with the addition theorem of the spherical harmonics 



(21 + l)P,(cos 6') = 4tt £ Y* m (6 A , <p A )Y lm (e , if) , (76) 

m 

one readily derives m = mf — after the angular <£> integration in eq. flTBl . 

For simplicity we shall work out the correlation integral only for the 
transition from a circular to a circular state. The result is 

jrii-rif 

a if (A) « —= c (n t ,n f )V:C:;icos9 A )e^f- n ^ 



X 



oo 

2 







drr j ni _ n/ (Ar)i? n/n/ _i(r)i? nini _i(r) , (77) 



c(ni,n f ) 



J^dcos6r:j_::( C os6)V:;J(cos6)V:;i 1 1 (cos6) , (78) 



where V™(x) is the normalized, associated Legendre polynomial. We have 
taken the leading term in I sum, I — |n$ — n/|. In the large n limit of the 
principal quantum number, 

. ,. „.. . .. (T( fn. - n r 4- 3V9\ 1/2 

(79) 



c(n, , n/ ) ~ .-^(^1)1/4 f r((n t -n / + 3)/2 \ 1/2 



which is of order unity. The remaining radial integral is neither small for the 
region of 

2V2~tt 1 . fn 

an n 3 

with 5n = rii — rif <C riij. For the transition around n « 20, 



A <— ^-Vl' (80) 



- llmeV( — f' 2 , (81) 

rr n 

20 

A < A ma:E , A max = 370e\/(— ) 3/2 . (82) 

n 

Thus, the correlation is large for a range of A ^> Ei — Et. 
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5.2 Field assisted pair decay 



Suppose that strong microwave is irradiated to a Rydberg state which 
then decays into a neutrino pair + another Rydberg state |/); 

Ri^R f + ViVj . (83) 

Electronic transition \i) —>■ \ f) under a strong EM field may be dealt with 
adopting quantum mechanical treatment of multiphoton processes [12] , or its 
extention. One needs a theoretical formalism in order to properly incorporate 
strong field effects. 

We consider a time dependent EM field in the Hamiltonian taking the 
radiation gauge, 

Hli t)=M£ + A. (84) 
m e 2m e 

The plane-wave microwave of a linear polarization is given by the vector 
potential of the form, 

— * 

E 

A(t) = — mn(wt) , (85) 
uj 

while a circularly polarized case is 

— » Eq 

A(t) = y= (sm(ujt)e x ± cos(ut)e y ) , (86) 

where e*j is the orthonormal unit vector along i— axis. We took z— axis as 
the direction of light propagation. In both cases of the polarization E Q is the 
rms amplitude. 

In discussions that follow it is important to distinguish whether the field 
is in the strong or the weak range of the strength. This may be characterized 
by interaction strength relative to its frequency; 

eEnp , i En ,GHz^ 

m e (jj^ V cm 1 oj 
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Depending on a combination of parameters, E /(u 2 n), this can be either in 
the strong or the weak field range. Even if the field strength is strong in this 
sense, it may be arranged that the field is not large enough to ionize Rydberg 
electron, eE < a/{n 2 aB) (attractive force from nucleus), which means 



We would like to treat the weak interaction process alone perturbatively, 
and solve interaction of atomic electrons with microwave as analytically as 
possible. As discussed in Appendix B, the transition amplitude of the neu- 
trino pair emission is given by 



(A f (oo)\p f )(p f \U A (oo ^^HwUAih , -oo^ipilA^-oo)) . (89) 



Here H\y is the weak interaction Hamiltonian of the neutrino pair emission 
of momenta 

States taken as initial and final ones in eq. (l89i) . U A (t , — oo)|v4j(— oo)) and 
U A (t , oo)|.Aj(oo)) are bound state solutions of the Schroedinger equatioin 
governed by the Hamiltonian of Coulomb potential plus the microwave field. 
Their momentum space representation U A (t , — oo)\pi) and U A (t,oo)\pf) are 
used here. The approximation, taken by Keldysh [TJ] and the one we shall 
also adopt, is to neglect the Coulomb interaction during the occurrence of 
weak process, and use the plane- wave solution under the periodic field, known 
as the Volkov solution ;T2l : the time dependent part in this approximation is 



E < 7 x 10 n n~ 2 Vcm 



-1 



(88) 




■oo 



— oo 



(p f \U A (oo ,t)\p f )(pi\U A (t , -oo)|pi) oc 




negelcting irrelevant phase factors. 
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The Keldysh approximation is valid if correction of the binding is small. 
This condition is worked out in [H], and in our case it leads to the frequency 
condition, as discussed in Appendix B, 

u < 2.1 x 10 6 n- 4 GHz, (90) 

with n the maximal principal quantum number during the transition. The 
field strength is also limited as in Appendix B. 

The weak Hamiltonian Hy/ of neutrino pair emission is translationally 
invariant, hence its matrix element contains the momentum conserving delta 
function; 

(p f \U A (oo ,t)H w U A (t , -oo)\pi) = 

( ( A. _ £_ _ El _ e 2) t + w zMi r~ gjjw 

\ 2m e 2m e m e 
x^|(27r) 3 5(pl - p f - q x - q 2 ) £<M^|0> ■ (f\j^) . (91) 

The neutrino pair is not observed, hence one takes the helicity and momen- 
tum summation of neutrinos. The helicity summation of 

EEKM^lo)-(/I^K)l 2 , 

hh' ij 

has been examined in Section III in detail, yielding in the leading approxi- 
mation of l/m e , the marix element squared of the form, (2 — m\vn2j{E\E^)\ 
eq. fTSS]) . times the electron wave function factors. 

To proceed further for computation of the transition amplitude squared, 
we insert a convenient identity, 

jd&J dE 12 5(A - q x - q 2 )5(E 12 - E 1 - E 2 ) = 1 , 



with Ei = J qf + mf neutrino momenta, and use 



/ / dK I ^2(2tt) 3 5(A - q x - q 2 )S(E 12 -E x - E 2 )f{q x , q 2 )\ 2 
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dA J dE 12 J ^f^5(A - ft - q 2 )8(E 12 - E l - £7 2 )j/(gi , g 2 )| 2 . (92) 



Neutrino momentum integration here is 



1 



, (mi+m 2 ) 2 \ / (mi-ra 2 ) 2 ^ 



( 27r ) 2 ^ \ S12 / V S 12 
x £ 2 2 (1 + ?-) - -^(1 + ^) 2 - mi m 2 

\ s 12 * s 12 

--(1 - (mi + m2)2 )(l - (mi ~ m2)2 )) , s 12 = E\ 2 - A 2 . (94) 

6 S12 S12 y 

In the massless neutrino limit, 

^(^i 2 ,A)^^( J E 1 2 2 -^). (95) 
The threshold behavior of this quantity at s 12 — > (mi + m 2 ) 2 is 

The transition probability is further simplified first by using 

(p1|e-^/( 2 ^)|A(-oo)) = (pHe-^IAl-oo)) = A(-oo)> , (97) 

and a similar relation for the final state, to replace the time dependent factor 
to e - l ( E i- E f) t , When the Fourier transformation back to the configuration 
space is made, one obtains 

/ d ~i0r m ~ pf ~ *H A f(o°mm\M-oo)) 
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The transition probability is then 

\(S-l)M 2 = ^f / dA\a tf (A)f dE 12 K*f(E l2 ,A) 

dtF lf {E 12 ,A;t)\ 2 , (99) 

T if (E 12 ,A;t)=exp[-t\A lf t+f t dh eA ' A{tl) )], (100) 

\ J —oo 17l e I 

A if = E i -E f -E 12 . (101) 

In the formula ( |99l) three important factors are separated; the inegrated 
neutrino factor K-j , the initial and the final electron factor (A) | 2 , and the 
time dependence factor Tif related to microwave irradiation. This essential 
simplification owes to the Keldysh approximation. 

The time integral for this S-matrix element invloves the time interval of 
infinite duration. In practice, it is important to understand a finite time 
integral of the form, 

f dhFifiEu^ih) = f dh exp[-z 1^+ T dt 2 &A ' A ^ ) 1102) 

J tg J to \ J —oo 7Tl e I 

The analysis of this time dependent phenomena is separated into two 
parts; in the first part we present the conventional multiphoton picture, 
which corresponds to the narrow band region of a more general analysis in 
the last subsection. There exists the additional wide band region which typi- 
cally exhibits the exponentially growing instability, a phenomenon very time 
dependent. Only a lower limit of the growing rate is estimated in the wide 
band region. The real process goes with these two mechanisms entangled, 
hence is complicated. 



5.3 Multiphoton picture 



31 



We postpone a general analysis of eq. (ll02p . and concentrate here on an 
approximate expansion that leads to interpretation based on multiphoton 
processes. The multiphoton process has been described in many textbooks; 
for instance in [12] for strong laser field, in [1] and (15] for microwave or rf 
fields. 

We shall first discuss the case of linear microwave polarization. Using the 
expansion in terms of the Bessel function Jn{%), 

e -ifesin^ = y jf^-iNut h= e J±l_^_ ( 103 ) 

eq. (11021) becomes a sum of simple exponentials. The time integral is then 
readily computed, leading in the large time limit to 

/ ^1^/(^12,^,^/;*!) -> 27iJ2JN(b)5(A if + Nuj). (104) 
J to N 

The argument of the delta function implies the simple energy conservation, 
due to ffTOll . 

It would be instructive to numerically estimate the important quantity 
that appears in these formulas; the magnitude of microwave interaction, 

6 ~^«0.2-^-^(™!)>. (105) 

The momentum scale A has been set here to around 2 times the neutrino 
mass. This parameter h can be very large. The Bessel function Jat(6) is 
maximal at b ~ \N\. 

Using the standard formula, 

MmM (A)P 

t->cc t V 1 

we derive a time independent rate w = lim^oo P(t)/t of the form, 

eA • E 



m P uj 2 



m = J2 w N, w N = irG 2 F dA\a if (A)\ 2 \J N 

N J 

roc 

/ „ dE 12 5(E 12 — Nuj — Ei + E f )K*f(E 12 ,A) . (106) 

J ^ A 2 +(m 1 +m 2 ) 2 
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There is a minimum number of photons N in the N summation; N > N 
with 



This N can be negative. A negative N means that stimulated microwave 
photons are emitted (instead of absorbed) to cause the neutrino pair emission. 
The requirement for a negative N is 



For the weak field the neutrino pair emission rate is in proportion to 



Interpretation of this result is that there are contributions from N photon 
absorption for N > or — N emission for N < 0, which feed or take away 
energy \N\u to cause the multiphoton transition. Roughly, the relation to 
the neutrino mass mi+m 2 ~ Ei — Ef + Nu holds. Suppose that the final 
state is specified. When Ef > Ei (the case of upper level), only microwave 
absorption is possible for the pair emission. When Ef < E { (the case of 
lower level), both absorption and emission is possible according to the sign 
of mi + m 2 — (Ei — Ef). Since the rate is maximal at \N\ ~ b oc E , the 
adjustment of the field amplitude Eq can help to locate the position of the 
threshold of Ei — Ef = m x + rri2- 

The neutrino pair emission accompanying iV(> 0) microwave absorption 
occurs as if a hypothetical heavy 'boson' of mass Nuj + Ei — Ef decays 
according to the rate wn- Let us clarify this process in the weak field limit. 
The neutrino pair emission caused by N multiphoton transition occurs with 
interaction strength, 




(107) 



Ei — Ef > mi + m 2 ■ 



(108) 




(109) 



G F 1 eA-E { 



\ N 



(110) 
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This is the El transition repeated N times plus the weak pair emission. The 
hypothetical 'boson' does not have a definite momentum to be transimitted 
to the neutrino pair, but the conservation law is replaced by 

(27r) 3 5(A)^MA)| 2 . 

There is nothing special about this, because both initial and final states 
are not momentum eigenstates. The momentum distribution has a width 
w a 2 m e /n 3 . The total rate is a sum over N of many multiphotons. A large 
mass 'particle' of iV 3> 1 might be called a heavy electron due to many 
photon clouds. The rate of a very heavy electron decay is suppressed by 

J_ 2 E i -E f -E 12 2N 

[ N\ } { 2lo ] ' 1 ] 

because the Bessel function behaves as |Jjv(6)| ~ (b/2) 2N /(N\) 2 . 

To obtain a large rate, it is necessary to have many contributions of 
different N. The property of the Bessel function J/v(&) tells that the N sum 
can be large when N is of order b or less. The ratio b/N, in particular 

b eA-E 



\ N o\ ^Ja 2 + (mi + m 2 ) 2 m e u ' 
can however be small, especially if A mi + m,2, because in this case 



112) 



' 0[^]~ 0,0.9x10-]-^^. (113) 



| iVo | m e uj Vera 1 uj 

Thus, it is important to have for N a small, or better, a negative number 
such that the iV sum contains contributions from small TV's. 

A better, but still crude way of estimating a number of multiphoton 
contributions of many paths, hence the pair emission rate is as follows. Take 
Eq along z— direction. For a given A z , there is a region of relatively large 
value of the Bessel function Jjv(&) at b ss iV 3> 1, 

^(^)~c(^)V3, c=i^~0.45. (114) 
v m e u 2 ' eA z E J 2 2 / 3 3 1 /6 7r v ; 
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Replacing the Bessel fuction by this gives the rate, 

w = nc 2 G 2 F f dAK^C-^ + E i -E f , A)|a 4/ (A))| 2 (^) 2 / 3 . (115) 

The threshold behavior of the rate based on ( 11151) may be derived using 
the neutrino factor Kf- near the threshold, (|96l) . For this estimate we assume 
the correlation integral of order unity ay = 0[1], to give 

(ma + m 2 ) 2 eE ^ V 



(116) 



where J is 



J= /^(l -p 2 ) 5/6 [ dzz- 2/3 VT^, c 2 J~0.15. (117) 

JO JO 

The N summation is limited by 

N < N max , iV max = 0[ if\ ~ 0.2— - — — ( ) , (118) 

where Aj/ is a typical momentum transfer of order mi + m 2 . 

The enhancement factor R may be defined relative to the standard rate 
near the threshold, G F (mi + m2) 5 /(157r 3 ); 



_ 45vr 3 c 2 J 



X 



^ m e u 2 5/3 (m 1 m 2 f /2 
16 eE c<j 2 Ajj(m 1 + m 2 ) 7 

+ ^ - E f ) 2 - K + m 2 ) 2 ) . (119) 



For the weak field of eE /(m e uj) <C 1, the threshold appears at Ef = Ei — 
nil ~ m 2- It may also appear as a threshold of the field amplitude E . A 
large power 8/3 implies that the rate quickly increases once the threshold is 
passed. 
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Well above the threshold one may use eq. (195|) . and the rate and the 
enhancement factor become of order, 

27c 2 G|o; 5 m e u 2/3 E t -E f im 

448 [ eE } UW+ lu ' ' l ^ Uj 

p-s.e^) 14 ^^^ 2 -) 1 / 3 , (121) 

m e uj uj 

taking \f = mi + rri2- This factor is very sensitive to the microwave power 
oc P 7 / 3 and its frequency oc co~ 5 , and its precise determination requires a 
more elaborate computation. 

There is also contribution from N < b region of the Bessel function 
Jat(&) ~ J%. This contribution is estimated as 

WN<b ~ ° [ 240 7e^ ] ■ (122) 

We shall finally consider the case of circular polarization. In this case 

E A 

A ■ A{t) = —^z — sin #a sin(u;£ ± (p&) , (123) 

— * 

@a , ^Pa being angle factors of the momentum A. The expansion in terms of 
the Bessel function becomes 

exp[ze - W ] = E^ d ( kl 2 exp[-2iV>t T iN d <p A )] . (124) 



m P ur 



The pair emission rate for a circular to a circular transition is given in terms 
of the correlation integral by 

w = 2tt 2 G 2 f / dAA 2 / dcos0 A |M A )l 2 PnK-i(cos0 A ) x 
JO J—i 1 



\2 



X 



9 \ Ei - E f T (rii — n f)u) - \J A 2 + (mi + m 2 

\Jn t -n f ( ^2^ ) |2 ^ ( ^ T (W< ~ ^ ~ Ef ' A) • (125) 
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A detailed and more precise rate computation in the multiphoton picture 
shall be presented elsewhere, and be compared to a more general approach 
presented in the following subsection. The method described in the present 
subsection has a limitation, and is interpreted as a part of more general 
approach we shall now describe. 



5.4 Relevance of parametric resonance 

Consider the time dependent part of the rate, 

W if {t; E l2 , A) = Re (j? f {t) £ dt x ^/(ti)) , (126) 

which appears in the transition rate given by the time derivative of the tran- 
sition probability (1991) . namely 

w(t) = / dA\a lf (A)\ 2 / dE l2 Kff{E l2 , A)W if (t ; E 12 , A) . 

(127) 

We define a complex function G(t;a,b) in terms of new dimensionless vari- 
ables a , b 

G{t;a,b) = e ia ^+ ifesinw ' f dt' e -^t' -ibsm^t' ^ ^ 28 ) 

(129) 



o 

Ej-Ef- E 12 eA-E 

' 9 

on m P u; z 



such that Wi f {t ; E X2 , A) = G(t ; a(E 12 ) , 6(A) ). 

The multiphoton picture in the preceeding subsection corresponds to an 
infinite time limit ignoring a coherence in the computation of G(t;a,b), 
which leads to a time independent, constant rate w(oo). However, there 
exists an intrinsic instability in some parameter region of (a,b), which we 
now discuss. A coherence effect at finite times is crucial in this discussion. 
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The quantity (I129p satisfies a coupled differential equation, 

/ ^ + oj 2 (a + b cos ut) 2 -buj 2 smut \ ( ReG(t) 

\ buj 2 smLut + tu 2 (a + bcostot) 2 J y ImG(i) 

This system has an intrinsic frequency scale \a\u — \Ei — Ef — E\^\ such 
as given by the energy difference between the initial and the final states. 
This intrisic scale is further modulated by a periodic variation of parameters 
whose frequency is uj and amplitude is |6|. Cooperative effects of external 
modulation with the intrinsic property gives rise to interesting phenomena 
of the parametric amplification. 

In our time dependent problem, the initial condition is specified as 

G(0) = 0, G(0) = 1 (131) 

Thus, there is no way to avoid the instability of the parametric resonance, 
once the parameters (b , a) fall in the instability band. 

General theory [T7] of linear differential equations with periodic coeffi- 
cients indicates solutions of the Mathieu type, and the unstable and stable 
band structure appears in the parameter (b , a) plane. In the unstable band 
the exponential growth is observed; 

Re^^-^-Zf-^fiut)) , (132) 

with fir) a periodic function, hence the instability greatly expediates de- 
pletion of the prepared state. The parameter b is essentially oc the total 
momentum of the neutrino pair A z projected onto the microwave electric 
field direction taken z— axis here, which is also a typical momentum trans- 
fer in the electron transition \i) — > |/). Another one a cx E12 — Ei + Ef 
is the total neutrino energy minus the mass difference Ei — Ef. Thus, 
(6 , a) oc (A 2 , E 12 — Ei + Ef), and the unstable band structure in (b , a) plane 
signifies where in the phase space of the neutrino pair contributes to the 



0. (130) 
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emission rate. The instability signifies an exponential decay of the initially 
prepared state. 

The relative weight of b term in eq. (11301) . the magnitude signifies 
the importance of the parametric amplification. Roughly, the narrow band 
region is in the parameter region of \b\ <C \a\, and the wide band region is 
in |6| ^> \a\. The multiphoton picture explained in the previous subsection 
corresponds to a narrow band of instability. The iV— th band in the narrow 
band region corresponds to a mass difference of initial and final states, Nu + 
Ei — Ef shifted by the energy input of N microwave photons. The narrowness 
implies weaker rates. The diagramatic interpretation of the narrow band 
decay has been given in the literature [18J. The correspondance between the 
two regions is given by 

Y,\ J N(^ 1 ^)\ 2 ^(Noo + E i -E f -E 12 ) <-> W if (t;E 12 ,A). (133) 

Rather than a discrete N sum of multiphotons there is a continuous spectrum 
of heavy electrons of mass Eyi present in the wide band region. 

The wide band region gives a much more enhanced time dependent rate 
than the narrow band multiphoton result presented in the previous subsec- 
tion. If (b,a) lies deeper in the instability band, namely, the larger |6|/|a| 
is, the greater the rate is. The rate readily exceeds order unity (in the unit 
of l/ou, that is T — 0[u]) deep in an instability band. In the wide band re- 
gion there is no definite number of muliphotons N, and the continuum broad 
mass 'E[ — Ef range all contribute to the instability. If one experimentally 
arranges that there is no radiatively decaying state as in the cavity QED, 
then this means that the neutrino pair emission is expediated; enhanced pair 
emission. It is thus important to depict the structure of stability-instability 
bands in (b , a) plane. 

The phase space region in terms of (£12 , A) of the wide band region 
is estimated as follows. One can imagine that the most important region is 
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restricted to \a\ < \b\ due to an experience in the Mathieu equation. For large 
6's, the coupled equation (11301) approximately decouples, and the coefficient 
term of the form, oc (1 — cos2co>t) implies that the relevant (b , a) region is 
indeed deep in wide band regions. The exponential growth rate given by 
e Xu)t /(\uj) is of order A ~ 0.15/2, as shown in [23J. 

The two constraints on the wide band region \a\ < \b\ and the mass 
restriction correspond to a region of (£12 , A), 

(E 12 -Ei + E f f < (—) 2 A 2 Z , E\ 2 > A 2 + (m, + m 2 f . (134) 
m e u) 

This gives a neutrino pair momentum integration of order, 

2 f dA\a lf (A)\ 2 [ E ^ Ef+e{Azl dE 12 K^(E 12 ,A)W if (t;E 12 ,A) ,(135) 

J JEi-E f -e\A z \ 

for e < 1 with 

e=| — |- (136) 

m e u) 

For a stronger field, the phase space area is quite different. In particular, 
towards and above the critical strength E c , 

£ c = ^~l.lx 10 5 Vcm~ 1 -^- , (137) 
e GHz 

the momentum space integration is changed to 

r E t -E f +e\A z \ 



2 / dA\a lf (A)\ 2 [6{A t -A z ) / dE 12 

J \ J Ei-E f -e\A z \ 



rE l -Er+e\A z \ \ _ _ 

+9(A Z -A*) I dE 12 )K^(E 12 ,A)W lf (t;E 12 ,A). (138) 

J v /A 2 +(mi+m 2 ) 2 / 

In the 6 > 1 limit the second term is dominant, and 

A. « l - fa -E f + ^/A2+A2 + (m 1 +m 2 ) 2 ) . 

We shall make a crude estimate for the pair emission rate in the wide 
band region by making two assumptions. We first introduce an average rate 
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factor (Wif), and next compute the rate far away from the threshold region 
for which one may use eq. (1951) . The Eyi integration then gives 

3 

« JdA\a lf (A)f^- 2 \A z f(W lf ). (139) 
It is thus expected to obtain the rate larger than 

- k %§rO s < w "> (l40) 

~ 6 x 10-.- (O^WID) i^f^^f , (14 1) 



using A maa; of order, eq.([82]). 

Once the large rate is confirmed, one may go to the threshold region, in 
which one replaces the energy integral by 

r Ei-E f +e\A z \ 

/ , dE l2 K^(E l2 , A) 

(mim2) 3//2 



(£. _ ^ + e |A z |)V(^ - E f + e| A,|) 2 - A 2 - (m x + m 2 ) 2 . 

(142) 



This gives the threshold rate, with e = eE / '(m e u), 
G 2 F Al ax ( mi m 2 yf> {{eAmax _ + ma)2 . 5/2 ((w )w/10) _ 



67r (mi + m 2 ) 2 e 2 



Above a field threshold of (mi + m 2 = 2m,, 



(143) 



2m e um u _ x m u u 400eV 

SVcm — , (144) 



eA max 50meV GHz A 

the rate quickly rises to a rate of order 



w 



2407T m e uj' 



{W if ) , (145) 
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taking = m v . This rate is m u /A max times the rate much above the 
threshold. 

A large factor by (A ma:r /m„) 5 in these rates is due to a larger momentum 
spread in sharply locarized circular Rydberg states. 

The critical field strength for the circular polarization is a factor \J~2 
larger, but the phase space is also different. 

It is important to keep in mind that the coherence should be maintained 
during the microwave irradiation. If this is possible for a long time, one may 
expect a huge growth factor (Wj/). The ultimate bound on (Wj/) is derived 
by the unitarity argument in the following way [21]. The unitarity for time 
dependent process requires 



With a given time dependent rate Wfi(t), this is essentially a bound on the 
allowed time T. It is appropriate to parametrize our weak process with 



due to the energy dependence of the weak transition rii — > rif by the neutrino 
pair emission. The most important dependence is on the final state rif, 
as indicated. We ignored less important dependence of A , A on rif. The 
requirement of unitarity eq. (I146j) is then roughly 



For estimate of riQ we take the allowed lowest state for the pair emission, 
Ei- E f > 2m u . 




(146) 




(147) 




(148) 




< 2m, 



V 



(149) 
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Taking as an order of magnitude estimate no = ya 2 m e /(4m^) gives 

A e A- < 5 _lf „ 3 x ^l^i^^f , (150) 
auj 2nj z rii m v 

corresponding to the maximal rate, 

Ae" ( (^) 2 - lV < 3 x 10V 1 -^(^^) 5 ( (^f - lV , (151) 
\ nj J GHz m u \ n f J 

with A = 10. 

The actual rate would be very large, much larger than the number given 
in eq. (!14ip . but presumably less than the unitarity bound f)15ip . A practical 
rate might be limited by actual experimental conditions such as the loss rate 
of coherence of Rydberg atoms. Most relaxation time may be arranged to 
be much larger than 10/u; ~ 10~ 8 s(GHz/u) such that the real rate may be 
close to the unitarity bound. There appears a real possibility of measuring the 
neutrino pair emission process from circular Rydberg states if the background 
rejection is successful. 

It is beyond the scope of our present work to precisely locate and further 
exploit the wide band regions for the parametric amplification of the neutrino 
pair emission. Detailed study of this aspect will appear elsewhere. 



6 A few comments on experimental method 

Clearly, one needs a systematic study, both theoretical and experimental, 
to implement our idea of the neutrino pair emission from circular Rydberg 
atoms. We shall be content here to make a few, rather trivial, comments 
towards a more organized study. 
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Laser irradiated pair emission process is relatively straightforward, hence 
we shall focuss on the microwave irradiated Rydberg atoms. 

6.1 Key ideas for experimental success 

One has to avoid the danger of disapperance of the intial Rydberg states 
via ordinary radiative decay, since it might also be enhanced by the wide 
band parametric resonance. For this purpose, we may use circular Rydberg 
states as the initial prepared state, for which radiative decays are very much 
suppressed except n — ► n — 1 El transition. This El transition can be 
suppressed for instance by the cavity QED effect [2]; radiative decay mode 
within a cavity is modified by the boundary effect, and if the wavelength 
A > 2d, a size of the cavity, the decay may be inhibited. The inhibition 
prevails for all multiphoton transitions, since all photons in this case satisfy 
the same condition A > 2d once one photon transition is inhibited. One 
must however arrange experimental apparatus such that the cavity does not 
interfere the microwave irradiation |25j . 

Another important issue is how to unambiguously identify the pair emis- 
sion process. Identification of final states is most important for the determi- 
nation of the threshold, and this can be done by the field ionization technique 
PQ. The threshold for the neutrino pair ViVj appears at 

eAif + Ei - Ef = nii + rrij , (152) 

with e = eE / \m e u). 

The best way for unambiguous identification of the weak process is to 
measure parity violating (PV) effects which are absent in QED processes. 
The simplest of this kind is to use circularly polarized photon beam and to 
measure the difference of atomic transition rate between h = ±1 polarization. 
The asymmetry of the emitted photon distribution along the direction of 
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the polarization oc J ■ p y is another parity violating measurable. Detailed 
theoretical study of this effect will appear elsewhere. 

6.2 How to proceed 

Since the small rate is a critical issue, an organized strategy of experimen- 
tal efforts is important. We believe that the first step should be discovery of 
the neutrino pair emission from excited atoms, and then one should steadily 
approach the smaller energy scale towards the pair emission threshold. Pre- 
sumably, at a few times twice the heaviest neutrino mass of order 0.05eV, 
namely at energy ps 0.3eV, one can hope to observe a signature of the dif- 
ference between the Majorana and the Dirac neutrino. The final step is 
a precision neutrino mass spectroscopy along with measurement of mixing 
angles. 

The use of atoms for neutrino physics is a new concept. It is evident 
that many R and D are required for this project, but one can initiate both 
experimental and theoretical efforts by use of modest human and budgetary 
resources. This is perhaps the ideal way towards a difficult physics goal. 

1 would like to thank A. Fukumi, I. Nakano, H. Nanjo, and N. Sasao 
for many stimulating discussions, in particular on experimental feasibility of 
the ideas presented here. N. Sasao has also pointed out a few mistakes in 
the original version of this manuscript, which is much appreciated. I also 
appreciate Y. Okabayashi for providing Figure 2 in the revised version of 
this paper. 

Note added 

After submission of this paper for publication, we realized that laser ir- 
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radiated pair emission at lower thresholds 2m! and mi + m 2 is useful for 
detection of the relic cosmic neutrino of 1.9K via the Pauli blocking effect. 
See [27] for this observation. 

Note added in proof 

After submission of this paper for publication analysis based on the opti- 
cal Bloch equation has been performed, and the greatest enhancement factor 
when 2 laser is irradiated, adding another laser to |/*) «-> |/**) transition, 
has been obtained. This work will be published elsewhere. 



7 Appendix A; Basic formulas of Majorana 
field 

7.1 Dirac and Majorana fields 

Using the representation of the Clifford algebra 

7a 7/3 + 7/37a = ^9a(3 , (153) 

for 7 a (Greek a = , 1 , 2 , 3 and Roman * = 1,2,3), 

o)- 7i= U »)• (154) 

( — o i \ / -1 \ 

7o7i = ,76 = ^7o7i7273 = , (155) 

V o a J \ o l j 

*=| l), (^)=id-7 5 W, (°)=i(l + 7^, (156) 
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the Dirac equation 



h a d a iJ -mip = 0, (157) 

is written in the following 2-component form; 

(id t — ia ■ V)y? = mx , (id t + ia • V)% = rmp . (158) 

The identification, x — i^P*, in the Dirac equation gives the Majorana 
equation; 

(idt — io • V)<£> = ima^y* , (159) 

with m the neutrino mass. The 2-component spinor if belongs to (1/2 , 0) of 
the irreducible representation of the Lorentz group, while ia 2 (p* to (0,1/2. 

Explicit solution of helicity eigenstate is derived by solving the helicity 
eigen-equation of eigenvalue h — ±1 , 

^-»>(:)-°- (:)-*U + + «)- ^ 

The full plane-wave solution to the Majorana equation is then 

(p(x) = e~ ip - x NZ - e ip - x Ep ^ hp N*ia 2 C , (161) 

{ = «?,*)=( J + ). (162) 

y /i(pi + «p 2 ) / 

An equivalent form of solution is obtained by using 

E p + hp _ E p + hp 
m ^ E p — hp 

When the helicity operator — ia ■ V/|V| is applied, the first term gives the 
multiplicative factor h, while the second gives the factor —h. Thus, the 
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consistent quantum interpretation of particle annhihilation and anti-particle 
creation of two terms is given by eq. ([6]) . 

The free quantum Majorana field is described by the Lagrangian density 
01 of 

Cm = - ((y^<9 y? + (pUa ■ Vy?) + (h.c.) - im ^a 2 f* - (/9 T o- 2 ^)) • (163) 

The Majorana particle must be quantitzed according to the anti-commutation 
rule; 

{<p a (x,t),<pl(y,t)}+ = 8 3 (x-y)5 aP . (164) 
Mode decomposition in terms of plane-waves is given by 



<p(x) = 

p,h 



where 



c(p, h)u(p, h)e~ ip ' x + c f (p, -h), 



E p + hp 
\E p -hp 



-ia 2 )u*{p, hy px 



{c(p, h) , c\p , h')}+ = 6p#6hh' , 
{c(p,h) } ctf,ti)} + = 0. 

Relation of discrete and continuous momenta is given by 

( 27r ) 3 r(3)/^ ^ V- x „/ d3 P 



(165) 

(166) 
(167) 



6: 



P,P 



V 



S {3) (p-ff), T,fp = vJ^- 3 f(p), (168) 



with V the volume of the normalization box. 

The normalization of 2-spinor consistent with canonical anti-commutation 
eq. (11641) is derived as follows. Computing the anti-cummutator with the 
plane-wave mode decomposition gives 

{<p a (S ,t) ,(fl(y ,t)}+ = 
J2 [eW*-V\N(p, h)\\p + hp 3 )(p + ha ■ p) af3 

p,h 



Ej v — tip 



(169) 
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We have used the relation, 

f (P, h )i ] (P , h) = (p + hp 3 )(p + ha-p). 
The correct anti-commutation relation thus requires 



(170) 



N(p,h) 



E p — hp 



2\\pE p (p + hp 3 ) 
Useful relations on the normalized wave function, 



(171) 



«(p. h) = - 



E p - hp I p + hp 3 
2\pE P (p + hp 3 ) I h(p l +ip 2 ) 



(172) 



arc 



u(p,h)ri(jr,h) = hi-%)(i+h?'''- 

4 h,p 
u f (p, h)u(p, h) = ~(1 - ^). 



P 



(173) 
(174) 



In quantum field theory it is important to identify Hamiltonian, momen- 
tum, and propagator. They are given for the Majorana field, 



H(x) = -icp'a ■ Vy? + (h.c.) + — ((p^a 2 cp* - ^ T cr 2 f) 
H = J d 3 xH(x) = J2 E p J(p, h) c(p, h) , 

p,h 

V = J <i 3 x<^(x)(— iV)ip{x) = ^2 pc\p ', h) c(p, h) , 



(0\T(ip(y)^(x)) \0) = -a-d 
(0\T(<p{y)<p{x))\0)=imcr 2 



d A p e 



ip-(y-x) 



(27r) 4 p 2 — m 2 + ie ' 
d 4 p e ip < y - x) 
(27r) 4 p 2 — m 2 + ie 



(175) 
(176) 

(177) 

(178) 
(179) 



The energy ( 11 76 ft and the momentum fl 177ft formulas of the Majorana field 
establish the corecctness of identification of c'(p,h) and c(p,h) as particle 
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creation and annihilation operators. The second form of the propagator (11791) 
is characteristic of Majorana field that does not conserve the fermion number, 
and does vanish for the Dirac field. 



8 Appendix B; Weak perturbative process for 
Rydberg atoms under microwave irradia- 
tion 

Consider a Hamiltonian system Hq perturbed by 2 types of generally time- 
dependent interaction, one of which with H + V A (t) is treated as solvable, 
and Vb as very weak; 

H = H + V A (t) + V B , H = ^- + V c (x), (180) 
(id t -H -V A (t))\A n (t)) = 0, (181) 

\A n (t)) = U A (t,t )\A n (t )), U A (t,t ) =ex P H / d*i(iZo + VA(*i))], 

(182) 

(A n (t)\ = (A n (t )\U A (t ,t) = (A n (t' )\U A (t' ,t). (183) 

In the interaction picture, the weak vertex is 

V B (t) = U A \t , t )V B U A (t , t ) . (184) 

The transition S-matrix S is computed using solutions \A n (t)) of (I18ip . from 

(AfitoWB^lMto)) = (A f {Q\U A \t l ,QV B U A {t ll U)\A l {t Q )). (185) 

By taking the time limit of t' — > oo and to ~~ ^ — oo, the transition matrix 
element is to lowest order of Vb, 

(S-l) fi ~-i [°° dttiAfioo^U^ihiOo^BUAih^oo^Aii-oo)). 

J — oo 
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(186) 



Expand states at time = ±00 in terms of a complete set of momentum 
eigenstates; 

]M _ x))= jfpjm^2)l. (187) 

(s - r )f , ~ -i£ *, / d% (A ff' m u A( «> ■ tOWfrft , -00) 

(2^ ' ( 188 ) 

We shall need to compute the transition matrix element for momentum eigen- 

st atesj 

(^1^(00,^)1^^1,-00)1^), (189) 
sandwitched between momentum state wave functions of 

(A f (oc)\p f ), (^lA(-oo)). (190) 

States here 

U A {t , -00) \pi) , (p f \U A (oo , t)) 1 " = U A {t , 00) |p» , (191) 

are solutions of the Schroedinger equation, eq. (j!8ip . 

These momentum state wave functions for H-atom are given in [20] • For 
instance, those of circular states of I = \m\ = n — 1 are 

^ n -Uv) - e^(l - cos 2 , (192) 

with m = ±(n — 1). For general states, using the atomic unit of po = 1/ob 
+„ta(p) = ^i(p)n m (9,^), (193) 

p -M«(d^ cK - ,( STT ) " (194) 
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where C^(x) is the Gegenbauer polynomial (of orfer N). 

Both states, {7a(£ , — oo)|Aj(— oo)) and ?7a(oo , £)IAf(°°))> contain the Coulomb 
and microwave field effects. Keldysh [13] neglects the Coulomb interaction. 
According to Reiss [H], this is justified if 

r = x— , 195 

m e cj 4 

is small. is the number of modes in the field). This reduces to 

r<-|^, (196) 
in our problem. E c is the maximum field allowd below ionization, given by 



a 

(n 2 a B ) 

Thus, the above parameter 

r < (a 2 ^) 2 \ ~ 4.3 x 10 5 (— ) 2 (H)8 . (i 98 ) 
uj n s uj n 

This parameter is small if 



uj < 2.1 x 10 6 GHzn-\ (199) 

and 

E <E C = a h/2 m 2 e n' A ~ 5.2 x lOVcm" 1 . (200) 
Thus, we arrive at 



U A (t,-oo)\d~eM-i t dhi^ + ^^ + ^T^-m, (201) 

J-oo) 2m m 2m 



d) 2m 

which is the Volkov solution [12], the exact plane- wave solution under a 
periodic field. Using an explicit form of A(t), 

E n 

lt ) 
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/ dtiAitt) = — mn(wt) + (const.) . (202) 
J t ) UJ 



Using the Volkov state, (neglecting irrelevant constants), the matrix ele- 
ment is 

(p f \U A (oo , t)V B U A {t , -oa)\pi) 

-expH f t + ® " &l " Eo sin(ujt) ) ] <g>|Vb|^> ■ (203) 

I 2m mo; 2 / 

Time dependence of the relevant matrix element of the integrand of H188[) 

(^/(^^(p/IC/^oo^JWaC^-^I^^IA-C-oo)), (204) 
is in eq. (12051). 

High frequency Floquet theory (HFFT) of [26J that predicts dressed 
Coulomb potential caused by averaged field irradiation appears irrelevant 
in our paremeter range of field strength and frequency. 
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